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Abstract 

We study the dynamics of the Forced Logistic Map in the cyhnder. We compute a 
bifurcation diagram in terms of the dynamics of the attracting set. Different properties of the 
attracting set are considered, as the Lyapunov exponent and, in the case of having a periodic 
invariant curve, its period and its reducibihty. This reveals that the parameter values for 
which the invariant curve doubles its period are contained in regions of the parameter space 
where the invariant curve is reducible. Then we present two additional studies to explain this 
fact. In first place we consider the images and the preimages of the critical set (the set where 
the derivative of the map w.r.t the non-periodic coordinate is equal to zero). Studying these 
sets we construct constrains in the parameter space for the reducibility of the invariant curve. 
In second place we consider the reducibility loss of the invariant curve as codimension one 
bifurcation and we study its interaction with the period doubling bifurcation. This reveals 
that, if the reducibility loss and the period doubling bifurcation curves meet, they do it in a 
tangent way. 
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1 Introduction 

We focus on the study of the quasi-periodically (q.p. for short) Forced Logistic Map (FLM for 
short). The FLM is a two parametric map in the cylinder where the dynamics in the periodic 
component is a rigid rotation and the dynamics in the other component is the logistic map plus 
a quasi-periodic forcing term. This map appears in the literature in different contexts, usually 
related with the destruction of invariant curves. For example, in [13] it was introduced as an 
example where SNAs (Strange Non-chaotic Attractors) were created through a collision between 
stable and unstable invariant curves. Since then, different routes for the destruction of invariant 
curves have been explored for this map, for instance see [25] and references therein. Some other 
recent studies are [H [21 [3 [lOl [E] [20] . 

On the other hand, the FLM is also related to the truncation of period doubling cascades. It is 
well known that the one dimensional logistic map exhibits an infinite cascade of period doubling 
bifurcations which leads to chaotic behavior. Moreover this infinite cascade extends to a wider 
class of unimodal maps. But when some q.p. forcing is added, the number of period doubling 
bifurcations of the invariant curves is finite. This phenomenon of finite period doubling cascade 
has been observed in different applied and theoretical contexts. In the applied context it has 
been observed in a truncation of the Navier-Stokes flow [111 [28] or in a periodically driven low 
order atmosphere model [6j. In the theoretical context, it has also been reported in different 
maps which were somehow built to have period doubling cascades [S] [21] , and more recently in 
the analysis of the Hopf-saddle-node bifurcation [7J . Actually, in |21j the FLM itself is given as 
a model for the truncation of the period doubling bifurcation cascade. 

The study presented here is more concerned with the mechanisms which cause the truncation 
of the period doubling bifurcation cascades than with the possible existence of SNAs for this 
family of maps. Concretely, we show (numerically) that the reducibility has the role of confining 
the period doubling bifurcation in closed regions of the parameter space. In the remainder of 
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this article we focus on the shape of this reducibihty regions. In |17l \TE[ ITU] we wih use the 
reducibihty loss bifurcation to study the self renormalizable properties of the bifurcation diagram 
and how the Feigenbaum-Collet-Tresser renormalization theory can be extended to understand 
it. See also |26j for a united exposition of the present paper with the other three cited before. 

This paper is structured as follows. In section [2] we review some concepts and results concerning 
the continuation of invariant curves for a quasi-periodic forced maps. We also look at the 
concrete case when the map is uncoupled. 

In section [3] we focus on the dynamics of the FLM. First we review some computations which 
can be found in the literature. Then, we do a study of the parameter space in terms of the 
dynamics of the attracting set of the map. For this study different properties of the attracting 
set are considered, as the value of the Lyapunov exponent and, in the case of having a periodic 
invariant curve, the period. Differently to other works, in our study the reducibihty of the 
invariant curves has been also taken into account. This reveals interesting information, for 
example we observe that the parameter values for which the invariant curve doubles its period 
is contained in regions of the parameter space where the invariant curve is reducible. The 
subsequent sections are developed with the aim of understanding the results presented in this 
section. 

In section [4] we consider the images and the preimages of the critical set (this set is the set where 
the derivative of the map w.r.t the non-periodic coordinate is equal to zero). We also consider 
the continuation in the parameter space of the invariant curve which comes from one of the 
fixed points of the logistic map. Doing a study of the preimages of the critical set we construct 
forbidden regions in the parameter space for the reducibihty of the invariant curve. In other 
words, we give some constrains on the reducibihty of the invariant curve. 

In section [5] the reducibihty loss of the invariant curve is considered as a codimension one 
bifurcation and then we study its interaction with the period doubling bifurcation. The study 
done here is not particular for the FLM. In this section we also give a general model for the 
reducibihty regions enclosing the period doubling bifurcation observed in the parameter space 
of section 13.21 



In section [6] we summarize the results obtained in the previous ones. We analyze again the 



bifurcations diagram obtained in section 3.2 taking into account the results obtained in sections 
landU 



2 Invariant curves in quasi-periodically forced systems 

In this section we briefly review some of the key definitions and results on the theory of invariant 
curves in quasi-periodically forced maps. These definitions will be useful for the forthcoming 
analysis of the dynamics of the FLM. 

2.1 Basic definitions 

A quasi periodically forced one dimensional map is a map of the form 

F : TxM TxM 

M ^ ( 

xj \fie,x) 
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where / E C'iT x M, M) with r > 1 and the parameter w G T \ Q. 

Given a quasi-periodically forced map as above, we have that it determines a dynamical system 
in the cyhnder, exphcitly defined as 



x = f{e,x) 



(2) 



Definition 2.1. Given a continuous function u : T — t- M we will say that u is an invariant 
curve of ^ if, and only if, 

u{9 + io) = f{e,u{0)), V^gT. (3) 
The value cj is known as the rotation number of u. 

An equivalent way to define invariant curve, is to require the set {{0,x) G T x M| x = u{6)} to 
be invariant by F, where F is the function defined by ([T]). 

On the other hand, if we consider the map F'^ we have that it is also a quasi-periodically forced 
map. Given a function u : T — )• M, we will say that « is a n-periodic invariant curve of F if 
the set {{9,x) G T x M| x = u{9)} is invariant by (and there is no smaller n satisfying such 
condition) . 

Since a periodic invariant curve of a map F is indeed an invariant curve of -F", any result for 
invariant curves can be extended to periodic invariant curves. 

Given x = uq{9) an invariant curve of ([2]), its linearized normal behavior is described by the 
following linear skew product: 

X = a{e)x, J ^ ' 

where a[9) = Dxf{9,UQ{9)) is also of class C"', x G M and ^ G T. We will assume that the 
invariant curve is not degenerate, in the sense that the function a{9) is not identically zero. 



Definition 2.2. The system Q is called reducible if, and only if, there exists a change of 
variable x = c{9)y, continuous with respect to 9, such that Q becomes 

y = by, } 

where b does not depend on 9. The constant b is called the multiplier of the reduced system. 

In the case that a(-) is a C°° function and co is Diophantine (see Proposition 1 in \20\). the skew 
product Q is reducible if, and only if, a(-) has no zeros, see Corollary 1 of [20]. Actually, the 
reducibility loss can be characterized as a codimension one bifurcation. 

Definition 2.3. Let us consider a one-parametric family of linear skew-products 

9 = 9 + u, I 
X = a{9,fi)x, J 

where uj is Diophantine and fx belongs to an open set of M and a is a C°° function of 9 and ^u. 
We will say that the system ^ undergoes a reducibility loss bifurcation at /Uq if 
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1. a{-,fi) has no zeros for /i < /io, 



2. a{-,fi) has a double zero at 6*0 for // = //q, 



3. J^a(eo,^o)/0. 



On the other hand, consider a system hke ^ with / a function, which depends (smoothly) 
on a one dimensional parameter (/ = /^). Assume also that we have an invariant curve u = 
of the system. We will say that the invariant curve undergoes a reducibility loss bifurcation if 
the system Q associated to the invariant curve {a{9) = a{9,fi) = Dxf^{6,u^{9))) undergoes a 
reducibility loss bifurcation as a system of linear skew-products. 

Given a map like ^ we have that, due to the rigid rotation in the periodic component, one of 
Lyapunov exponents is equal to zero (see |2]). Then the definition of the Lyapunov exponent 
can be suited to the case of linear skew-products as follows. 

Definition 2.4. If G T, we define the Lyapunov exponent of (|4| at as 

(7) 

We also define the Lyapunov exponent of the skew^ product Q as 



X{9) = limsup — In 

n— >oo 1^ 



n-1 

lla{9 + ju] 

j=0 



A = [ ln\a{9)\d9. 
Jo 



(8) 



If A is finite then, applying the Birkhoff Ergodic Theorem we have that the limsup in ([7| is in 
fact a limit and X{9) = A for Lebesgue a.e. G T. If a{9) never vanishes, the limsup in ([T]) is 
again a limit and coincides with A, but now for all 9 €T. 

Now, consider a map like ([T]) . If there exists an invariant compact cylinder where / is monotone 
and has a negative Schwarzian derivative with respect to x, then Jager has proved the existence 
of invariant curves, see [15] for details. On the other hand in [20J a result on the persistence of 
invariant curves is given, in terms of the reducibility and the Lyapunov exponent of the curve. 



2.2 Quasi-periodically forced maps which are uncoupled 

In this subsection we turn our attention to the maps which are of the same class of the FLM, 
in the sense that the quasi-periodic function can be written as a one dimensional function plus 
a quasi-periodic term. 

Definition 2.5. Given a map like Q we will say the the map is uncoupled if / does not 
depend on 9, i.e. f{9,x) = f{x) 

Note that if a map F is uncoupled, the also does. 

Proposition 2.6. Let he a one dimensional family of maps like such that for a fixed 
value Ho the map F^^ is uncoupled, that is F^^{9,x) = {9 + oj , f ^^^{x)) . Then any hyperbolic fixed 
point of ffj_Q extends to an invariant curve of the system for fi close to hq. 
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Proof. Suppose that there exists a hyperbohc fixed point xq G M of -F^p. We have that it can 
be seen as an invariant curve : T — )• M of F^^, with uo{9) = xq for any €z T. The skew 
product ([4]) associated to the invariant curve has as a multipher a{9) = ^f{9,xo) = f'fj,gixo), 
which actually does not depend on 9. Concretely we have that the system is reducible. 

Now we can apply the theory exposed in section 3.3 of [2D]. Assume first that /^„(xo) / 0. We 
have that a curve is persistent by perturbation if 1 does not belong to the spectrum of the transfer 
operator associated to the curve. Since the curve is reducible we have that the spectrum is a 
circle of modulus \a\. Using that the fixed point xq is hyperbolic we have that \a\ 7^ 1, therefore 
1 does not belongs to the spectrum. When /^^(a^o) = we have that the spectrum collapses to 
0. Then 1 does not belong to the spectrum of the transfer operator either. □ 

Note that, by considering for the periodic case, the result extends to any periodic point of 
the uncoupled system. 

This last proposition can be also proved using the normal hyperbolicity theory [14J. But this 
theory is only valid for diffeomorphisms, then the case of fL„{xo) = is not included. 



3 The Forced Logistic Map 

The FLM is a map in the cylinder T x M defined as 

= rUe) = 9 + u;, \ 

X = fa,e{^,x) = ax{l-x){l + ecos{2TT9)), J ^ ' 

where (a, e) are parameters and uj a fixed Diophantine number (typically in our study it will be 
the golden mean). 



3.1 Basic study of the dynamics 

Note that the FLM ^ is a q.p. forced system like Q, which depends on two parameters a and 
£. Moreover, we have that the function /q^^ which defines the map can be written as a logistic 
map plus a q.p. forcing term. In other words, we have that 

where ia{x) = ax{l — x) (the logistic map) and h£^a(9,x) = eax{l — x)cos(27r^) (which is zero 



when e is). Then proposition 2.6 is applicable to the map. In some cases it will be convenient to 
work in a compact domain. In this case note that when < a{l + |e|) < 4 the compact cylinder 
T X [0, 1] is invariant by the map. 

The one dimensional logistic map is known to exhibit a period doubling bifurcation cascade, 
that is a sequence of infinitely many period doubling bifurcations accumulating to a concrete 



parameter value. Applying proposition 2.6, we have that the period doubling cascade persists in 
the sense that, fixed a hyperbolic periodic orbit of period 2" of the logistic map, there exists a 
sufficiently small e such that the prescribed periodic orbit extends to a periodic invariant curve 
for e < eq. 

A natural question to ask is what happens in the converse sense. That is, consider the parameter 
£ fixed (equal to a small value) and let the parameter a increase. What happens to the period 
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Figure 1: The attractor of the FLM for different parameters of a when e is 0.01. The 
vertical axis corresponds to the x variable and the horizontal one to 9. The values of a (reading 
from left to right and top to bottom) are 2.4, 3.35, 3.5, 3.522, 3.527 and 3.529. The values of 
the Lyapunov exponent for each attractor are plot in figure [2j 



doubling bifurcation cascade of the logistic map? The answer to this question can be easily 
guessed with some simple computations. 

In figure [T] we have plotted the attractor of the FLM for several values of a when e is 1/100. 
We have also computed the Lyapunov exponent of the attractor, as a function of a in figure 
[2] The attractor has been computed by forward iteration of the map, using a transient and 
then plotting a certain number of iterates to obtain the approximation of the attractor. For the 
computation of the Lyapunov exponent we have used its definition as a limit (see ([T])). 

For the first three values of a we observe how the attracting invariant curve doubles its period 
twice. After two period doublings the curve begins to get more and more wrinkled until a 
Strang^ attractor seems to appear. This scheme is known as the truncation of the period 
doubling bifurcation and it is reproduced for any (arbitrarily small but fixed) value of e. 

In figure [2] we see that in the graph of the Lyapunov each period doubling of the attracting 
curve corresponds to a value of a where the exponent becomes zero. Note that between one 
period doubling and the next one the Lyapunov exponent of the system (as a function of a) has 
some critical values where its derivative goes to minus infinity. This behavior of the Lyapunov 
exponent is due to a loss and a recovery of reducibility of the attracting curve (see [20]). If one 
increases a, after a certain number of period doublings, there are no more bifurcations of this 
kind. Then, the attracting curve wrinkles until it becomes a strange set. In the graph of the 
Lyapunov exponent in figure [2] we observe that it goes to zero, but now it does it in a sharp 
way and then it becomes positive. This is a typical behavior also reported in other maps where 
there is a truncation of the period doubling sequence [21^ [28l [7| . Therefore the FLM has a lot 
interest as a simple model for this process. 

The destruction of the invariant curve shown above is indeed a subtle problem. For certain pa- 
rameter values the numerical computations of the attractor produce a set which seems strange. 

^ We will say that a set is strange if it is neither a difFerentiable curve nor a finite union of them. 
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Figure 2: Graph of the Lyapunov exponent (vertical axis) of the attracting set of the FLM ^ 
as a function of a (horizontal axis), for e = 0.01. The values corresponding to the attractors 
displayed in figure [T] are also shown. 



One might think that it is a SNA, but doing computations with higher accuracy, after several 
magnifications, the numerical approximation of the invariant curve does not seem strange any- 
more, but is an incredibly wrinkled invariant curve (see |121l20j). Actually, for certain parameter 
values one is not able to conclude (numerically) whether the curve is strange or not. 

In this last direction we have noticed that the detection of SNAs in the parameter space using 
the phase sensitive operator (for a description of the method see |23j ) can produce non-reliable 
results. In |20| it is given an example of a family of one parametric affine maps in the cylinder 
where the supremum norm of the derivative of the invariant curve goes to infinity when the 
parameter of the family tends to a certain critical value, but the invariant curve it still C°° 
before reaching the critical value. In this case, the phase sensitive indicator will fail and it will 
give a false positive for parameters close to the critical value. This kind of phenomena can 
happen when the method is applied to the FLM giving false positives of the indicator. 

To illustrate this statement we have reproduced the computation of candidates to be SNAs for 
the FLM done in In our computations we have considered several values of the parameter 
N (in the notation of j23]) when computing the phase sensitivity indicator. We have fixed a 
box in the parameter space of the FLM, then we have discretizated it in a rectangular grid 
of points. For each parameter in the grid we have computed the attractor of the map and its 
Lyapunov exponent. In the case of having negative Lyapunov exponent we have computed the 
phase sensitivity indicator. 

The results are shown in figure [3j The parameters values in white are the ones for which the 
iteration of the initial point diverged to —oo. The parameters in red correspond to positive 
Lyapunov exponent and the ones in blue correspond to negative exponent. The points in black 
correspond to the candidates to be SNAs obtained with the phase sensitive indicator. The 
different pictures correspond to different values of the order N. 

We can observe how the number of candidates decays when the order N of the method is 
increased. At the end we have much less candidates that what appeared to be in the original 
estimations of |24j . This is due to the fact that the indicator can not distinguish between 
a strange set from a very wrinkled (but smooth) curve. Then when the number of iterates 
increases, more and more candidates to SNAs are discarded. 



8 




2 2.2 24 26 28 3 3.2 34 36 2 22 24 26 28 3 32 34 36 




2 22 24 26 28 3 32 34 36 2 22 24 26 28 3 32 34 36 



Figure 3: Parameter space of the FLM The axes correspond to the parameters a (horizon- 
tal) and e (vertical). Candidates to SNAs are displayed (in black) for different orders (N) of 
computation, see the text for details. The values of taken for the computation (from left to 
right and top to bottom) are 10^, 10^, 10^ and 10^. 



Color 


Dynamics of the attractor 


Black 


Invariant curve with zero Lyapunov exponent 


Red 


Chaotic attractor 


Dark blue 


Non-chaotic non-reducible attractor 


Soft blue 


Non-chaotic reducible attractor 


White 


No attractor (divergence to — oo) 



Table 1 : Color coding for the figure |4j 



3.2 Parameter space and reducibility 

As in the previous section, we consider a certain subset of the parameter space {a,e) of the 
FLM. The parameters are classified depending on the dynamics of the attracting set, but now 
different properties of the attractor are considered. There are similar computations done in 
the literature in terms of the Lyapunov exponent of the attracting set and, in case of having 
a periodic invariant curve, its period [5]). In our analysis, when we have an invariant 

periodic curve we also take into account whether the invariant curve is reducible or not. Our 
computation of the zero Lyapunov exponent bifurcation curve has been done approximating the 
invariant curve by its truncated Fourier series, instead of approximating the rotation number uj 



9 



0.2 



0.15 



0.1 



0.05 




2 2.2 2.4 2.6 



2.8 



3.2 3.4 



3.6 



0.03 



0.025 



0.02 



0.015 



0.01 



0.005 




3.25 3.3 3.35 3.4 3.45 3.5 3.55 3.6 




3. 53.51 3. 523.533.543.553. 56 3. 573.! 



Figure 4: Diagram of the parameter space of the FLM (0) for a; = The axis correspond 

to the parameters a (horizontal) and e (vertical). For the correspondence of each color with the 
properties of the attractor see table [T] 



by rational approximations like in [22j 



3.2.1 The procedure 



We have considered three different rectangular subsets in the parameter space and we have 
discretized the subset in a grid of points. For each parameter in the grid we have computed the 
attracting set by forward iteration. Then we have computed the Lyapunov exponent of each 
attracting set. To compute this we have approximated the limit ([7]) for large values of n. We 
stopped when the variation between the estimation for n, |n and was lower than 10~^. 

For the parameter values with negative Lyapunov exponent we assumed that we have a periodic 
invariant curve and we have checked if it is reducible as follows. Given a initial value S TT, we 
can approximate the value of the invariant curve u at ^-s u{6o) ~ f^^{6o — Muj, xq) (for any 
xq). This allows us to compute the attracting curve in a mesh of points. On the other hand, 
we have that (if / is C°° and u Diophantine) an invariant curve is reducible if, and only if. 
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Dxf{-,u(-)) has no zeros (see Corollary 1 of [20j). Then we can use the discrete approximation 
of the curve to check this condition. 

The results are shown in figure |4j The parameter values where the attracting curve of the map 
has zero Lyapunov exponent are also displayed . In practice, this corresponds to the period 
doubling of the attracting invariant curve. The numerical computation of these bifurcations has 
been done as follows. 

Given u : T — t- M an invariant curve of ^ we have that it can be approximated numerically 
by its truncated Fourier series at order N. For more details on how this can be implemented 
see [SI [16]. Basically, we compute the Fourier coefficients as the zero of a suitable function 
G : M^iv+i _^ M2Ar+i_ ^.j^g particular case of the FLM we can add the parameters a and e 
as unknowns and extend to a function from R^^+^ to R^^+-'^. On the other hand, assume 
that we have an invariant curve which is reducible. Then we have that the Lyapunov exponent 
of the curve A(n) given by ([S]) is indeed differentiable (see |2Uj). If we have a curve u with zero 
Lyapunov exponent, we can use the function G = {G, A) as a continuation function for the zero 
Lyapunov bifurcation curve. Then the bifurcation curve in the parameter space is computed as 
the projection in the (q, e) coordinates of the set of zeros of G. 

To compute this bifurcation curve numerically it can be used a standard continuation method 
(see 127]). In our case the kind of continuation that we do is rather simple. Suppose that the 
bifurcation curve in the parameter space is regular (no critical points). Then at every point in 
the parameter space the curve can be expressed locally with one of the parameters as a function 
of the other. When one of the parameters is fixed, the other can be computed through a standard 
Newton method. To do the continuation of the bifurcation curve we can vary slightly the fixed 
parameter and compute the corresponding value of the curve. The selection of which parameter 
is fixed has been done depending on the estimated inclination of the curve in the parameter 
space. 



3.2.2 Description of the results 

In figure [4| we show the results of the computation and we have also plotted two successive 
magnifications of certain subsets of the parameter space. These regions have been marked with 
a dashed line in the picture. The points in the parameter space have been codified in different 
colors depending on the dynamics of the attractor, see table [Tj 

At this point let us recall the dynamics of the logistic map. We have that for certain values of 
a the map has a cascade of period doubling bifurcations. This is, we have that there exist a 
sequence of values with lim dn = d* , such that for a G {dk,dk+i) the logistic map 

n— >oo 

has an attracting periodic orbit of period 2^ and at the value a = dk+i the attractor doubles its 
period (from 2^ to 2^"^^). Let us recall that the logistic map is unimodal for any value of a (i.e. 
it has a unique point where the derivative of the map is equal to zero) . In the period doubling 
cascade, we also have that the attractor crosses the turning point between one period doubling 
and the next one. In other words, there exist values G (dfc,c^fc+i) for which the attracting 2^ 
periodic orbit of the map is the critical point. 

In figure |4] we can observe that from every parameter value (a,e) = {dk,0) it is born a curve 
where the attractor has zero Lyapunov exponent. Let us denote by Dk each of these curves in 
the parameter space, which corresponds to the period doubling from period 2^ to 2^^^. In the 
bifurcation diagram at the top of figure |4| there are plotted the curves Di and D2. In the bottom 
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left one there are displayed the curves D2 and D3, and finally in the bottom right one there are 
the curves D3 and 1)4. 

In the previous figure it has also been plotted the reducibility and the non-reducibility regions. 
We can also observe that from every parameter q a "cone of non-reducibility" is born, in the 
sense that there exist two curves in the parameter space which define a zone where the attracting 
invariant curve of the map is not reducible. Let us denote by C~ (and respectively by C^) the 
left (respectively right) boundary of the non-reducibility region born at the point a. 

We can observe how the curves C~_^ and define an enclosed reducibility region which contains 
the curve Di. Moreover these three curves seem to meet in a tangent way at the same point. 



3.2.3 Analysis of the bifurcation diagram 



Recall that in section 3.1 we have illustrated a truncation of the period doubling cascade for 
a fixed value of the coupling parameter e, see figure [l} The results described above on the 
parameter space of the map agree with the behavior reported there. Looking at figure [4] we can 
observe that each period doubling is confined inside a reducibility region. Moreover, when the 
period is increased each of these regions get closer to the line e = 0. Then, if one fixes the value 
of e (arbitrarily small) and let a grow one should expect a finite number of period doublings. 
Fixing e at a prescribed value eo > corresponds to fixing a line e = £q in the parameter space. 
If the enclosed regions of reducibility get closer to e = when the period grows, at some point 
the regions will be below the line e = Eq. On the other hand, the shape of the reducibility 
regions also explains why we observe a reducibility loss and afterwards a reducibility recovery 
between one period doubling and the next. At the same time, the bifurcation diagram displays 
many interesting phenomena which can be studied. 

The first phenomenon that can be observed in the bifurcation diagram is the birth of non- 



reducibility cones around the parameters values (a,e) = (cj,0). Recall that (definition 2.3) 
the reducibility loss of an invariant curve can be seen as a codimension one bifurcation. This 
codimension one condition defines a one dimensional curve in the parameter space, which is 
the boundary between the reducibility and the non-reducibility regions. With the notation 
introduced above, these boundaries will correspond to the curves and of the parameter 
space. 

Actually, the existence of a cone of non-reducibility around the points (a, e) = (cj, 0) is equivalent 
to prove that two reducibility loss bifurcations curves are born at these points. This will be 
proved under suitable conditions in [18], but the phenomenon can be explained heuristically as 
follows. 

Consider that we have an invariant (or periodic) curve u : T — )• R of the FLM. Using corollary 1 
of [20] we have that u is reducible if, and only if, Dxfa,e{'^ ^(0) no zeros. For the FLM, we 
have that the points {6, x) for which Dxfa,e{G, x) is equal to zero are those in the set {x = 1/2}. 
Therefore an invariant curve x = u{9) is reducible if, and only if, it does not intersect the set 
X = 1/2. Respectively, a periodic curve will be reducible if, and only if, none of its periodic 
components touches this set. 

On the other hand, when the system is uncoupled {e = 0) we have that the invariant (resp. 
periodic) curves are constant lines equal to the fixed (resp. periodic) points of the uncoupled 
system. For the logistic map we have that when the parameter a crosses the value cq (resp. q), 
the fixed points (resp. one of the components of the periodic orbit) crosses the set x = 1/2. 
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When one adds a small coupling to the system the invariant curve is no longer constant. Then 
when a is increased crossing the value cq (resp. Cj) the invariant curve (resp. periodic orbit) 
has to lose its reducibility, when it first touches x = 1/2, and then recover it again, when it is 
completely below (or above) x = 1/2. These two contacts of the invariant (or periodic) curve 
with the set x = 1/2 give place to the cone of non-reducibility when the parameters (e, a) 
are close to (0,co) (or (0, q)). Each boundary of this cone correspond to a loss of reducibility 
bifurcation curves, namely C~ and , in the parameter space. 

Another observable phenomenon in the bifurcation diagram of the figure |4] is that from every 
parameter value (a, e) = (di, 0) it is born a curve in the parameter space where the attracting set 
has zero Lyapunov exponent. Each of these curves correspond to a period doubling bifurcation 
of the attractor. For diffeomorphisms it is known (|5j) that the period doubling of an invariant 
torus is a codimension one bifurcation. Recall that the FLM is not invertible, therefore this 
theory is not directly applicable. Nevertheless, in the reducible case, one can apply a normal 
form procedure around the invariant curve, to obtain a one dimensional system in a neighborhood 
of the curve. We will formalize this argument in section [5] Then we can understand the period 
doubling bifurcation of the invariant curve as a period doubling of the reduced system (in the 
classical one dimensional sense). For the uncoupled case (e = 0) this happens at each parameter 
value a = di. Due to the codimension one character of the bifurcation a period doubling 
bifurcation curve is born at each of the parameters (a, e) = (dj, 0). 

One remarkable property (observed numerically) is that each period doubling bifurcation curve 
Di is confined in a reducibility region delimited by Sf_-^ and S~ . This can be explained analyzing 
the numerical procedure. Recall that the continuation of the bifurcation curve has been done as 
long as the estimated error was below a prescribed tolerance. On the other hand the function 
used for the numerical continuation contains the Lyapunov exponent of the invariant curve as 
one of its components. The method used to estimate the Lyapunov exponent behaves much 
worse when the reducibility is lost. Note that the Lyapunov exponent is obtained integrating 
numerically \a\Dxf{0,u{9))\. In the C°° case we have have that, in the reducible case the 
function is C°° whereas in the non-reducible case the function is not even bounded. This makes 
the numerical integration behave much worse when the reducibility is lost. Despite of that, we 
have the following reasons to believe that the invariant curve (at the parameter of bifurcation) 
is actually destroyed due to the reducibility loss. 

From the theoretical point of view, when the reducibility is lost the spectral problem associated 
to the continuation of the invariant curve changes drastically. When one tries to apply the IFT 
to an invariant curve there is a big difference from the reducible case to the non-reducible one 
([20]). Given an invariant curve and its skew product consider the transfer operator associated 
to the curve. We have that the invariant curve can be continued when 1 does not belongs to the 
spectrum of the transfer operator. In the reducible case the spectrum corresponds to the circle 
of radius equal to the exponential of the Lyapunov exponent of the curve. In the non-reducible 
case the spectrum corresponds to the disk of the same radius. When a period doubling occurs, 
there is an attracting invariant curve which becomes unstable, in other words its Lyapunov 
exponent crosses zero. In the reducible case, one might expect the spectrum of the transfer 
operator associated to the invariant curve (which is a circle) to cross 1 transversally when the 
bifurcation happens, the curve becomes unstable but survives. But in the non-reducible case 
there is no chance of transversal cross. 

From the numerical point of view, one can study the behavior of the invariant curve for 
the parameter values in the period doubling bifurcation just before losing the reducibility. 
In figure ^ we have displayed the invariant curve x{9) of the FLM ([9]) for the parameters 
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Figure 5: In the top left we have the invariant curve x{9) with zero Lyapunov exponent of the 

FLM ([9]) for the terminal point of the first period doubling bifurcation. In the top right and 

bottom left we have the first {x'(9)) and the second derivative {x"{9)) of the same invariant 

curve. In these pictures the horizontal axis corresponds to 9. In the bottom right we have (in 

the vertical axis on a logarithmic scale) sup (solid line) and sup |j;"(0)| (dashed line) as a 

6»eT eeT 
graph of the parameter e (horizontal axis), where is the invariant curve with zero Lyapunov 

exponent for each e. 



(a,e) = (3.3796,0.1423). This parameter corresponds to the last point in the invariant curve 
that we have obtained in the numerical computation. In the figure we also show the first and 
second derivatives of the curve. 

We can observe that the invariant curve x{9) seems to keep continuous and it does not fractaliz^ 
but what fractalizes is its derivative x'{9). Back to figure |4] we can observe that the period 
doubling bifurcation curve can be parameterized by the parameter e. Then for each e we can 
consider a{e) the parameter for the period doubling bifurcation curve and Xs(9) the invariant 
curve with zero Lyapunov exponent for the corresponding values of a and e. In figure [5] we have 
the supremum of |x£(0)| and |a;£(^)| as a graph of e. The figure indicates that sup |x"(^)| grows 

unbounded when we approximate certain critical value e ~ 0.1423. Moreover, looking at the 
graph of |x"(^)| this limit seems to be uniform when taking the supremum on any subinterval 
of T. 

This destruction process is similar to the fractalization process described in |20j . but the curve 
which fractalizes here is x'{9). The fact that the end point of the period doubling bifurcation 

^ Heuristically, we consider that a curve fractalizes if its length goes to infinity. 
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corresponds to the parameter where the invariant curve touches the Hne x = 1/2 was already 
observed in [22] . The computations of the bifurcation curve done in the cited work are based on 
the rational approximation of the rotation number u. Our method is based on the computation 
of the Fourier series of the invariant curve and it has the advantage that it allow us to compute 
the derivatives of the curves easily. 

To finish this section let us remark that this analysis is far from complete. In this direction, 
in the forthcoming sections |4] and [5] we present two studies in order to understand better the 
bifurcation diagram of the figure |4j In section [6] we summarize the different results and we 
analyze their implications on the cited bifurcation diagram. 



4 Obstruction to reducibility 



Recall that the set x = is always an invariant curve of the FLM. Apart from this set, when 
e = and a > 1 we have also the invariant curve given as Xa,o{d) = 1 — - . Using proposition 



2.6 



we have that for any a 7^ 3 this curve persists for small values of e. Let us denote by Xa,e the 
continuation of this curve. In this section we are concerned on the reducibility of this invariant 
curve. Considering the images and preimages of the set of points where the derivative of the 
maps is equal to zero, we will construct regions in the parameter space where the curve Xa,e can 
not be reducible. For instance we will see that a > 2 then |e| < 1 — ^ is a necessary condition 
for the reducibility of Xa,e- 

In this section we will consider different sets in the cylinder T x [0, 1]. These sets will be the 
closed graph of a curve or a subset of a graph. In general when we say that one of these sets is 
above (respectively below) of another, we mean that, for each value of G T the corresponding 
x-coordinate of the first set is bigger (resp. smaller) than the x-coordinate of the other set (for 
the same 6). The proofs in this section have been omitted because of their simplicity. 



4.1 A first constraint on the reducibility 

Definition 4.1. Given a q.p. forced map like ^, we define its critical set as the set of points 
on its domain where the derivative of the map (with respect to x) is zero. In other words, 

Po = {{e,x)eTxR\aj{e,x) = o}. 



In the case of the logistic map we have that Pq = {{9, x) G T x M] x = 1/2}. 

When we have a q.p. forced map like ([T]) which is C°° and oj is Diophantine, Corollary 1 of |20j 
implies that an invariant curve yo{0) is reducible if, and only if, 

Po n Graph(yo) = 0, 
where Graph(yo) = {iO,x) E T x [0, 1]| x = yo{0)}. 

In the particular case of the FLM when a(l + |e|) < 4 we have that the compact cylinder 
T X [0, 1] is invariant by the map. Then, any reducible invariant curve is either above or below 
the critical set. Actually, when the map is uncoupled {e = 0) we have that the invariant curve 
xo{0) = 1 — 1/a is above the critical set when a > 2. 
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Consider now Pi the image of the critical set by the map, namely post-critical set, which is 
defined as 



Pi = {{e,x) e T X [0, l]|6' = 61 + a;, x = f{e,x), for some (6',x) e Pq}. 

In the particular case of the FLM we have. 

Pi = G T X [0, l]|x = ^(1 + ecos(27r(^-a;)))| . 

Proposition 4.2. In the case of the FLM we have the following properties on the post-critical 
set when \e\ < 1. 

1. The set Pi is above the image of any other point (6, x) £T x [0, 1]. 

2. The points below Pi have two preimages in T x the points in Pi have one preimage and 
the points above Pi have no preimage. 

Consider the invariant cmve Xa,e of the FLM introduced at the beginning of this section. When 
e = and a > 2 we have that the invariant curve exists and it is above the critical set Pq. 



On the other hand, Xa^e is always below Pi due to proposition 4.2 Concretely, we have that 
a necessary condition for the reducibility of Xa,e is that Pq and Pi do not intersect. Otherwise 
there is no room for the invariant curve to exist without losing its reducibility. 

Then for any a > 2 it is necessary that, 

1 Q 

- < -(1 + ecos(27r(9)) for any 6* G T, 
which will be satisfied if, and only if, 

|e|<l-^. (10) 

a 

This gives a first constraint in the parameter space for the reducibility of the curve Xa,e- In 
figure [S] we have displayed this constraint together with the bifurcation diagram computed in 
section |3.2[ We can see that the constraint is quite sharp for values of a close to 2, but it 
becomes rapidly pessimistic for values of a far from 2. 



4.2 Further constraints on the region of reducibility 

In this section, we consider the preimages of the critical set to give additional restrictions to the 
region of the parameter space where the invariant curve Xa,e can be reducible. 

Definition 4.3. Given F a q.p. forced map like ^ we define P-k the k-th. pre-critical set as 
the k-th preimage of the critical set Pq. In other words, 

P_fe = {{e, x) G T X M| P^(e, x) G Co}. 

In figure [6] we have (in a solid line) the first four pre-critical sets for the FLM for a concrete values 
of the parameters. In the figure we have (in a dashed line) the post-critical set Pi. In the top 
left part of the figure we have displayed the sets Pq and P_i. Note that Pq lays completely below 



the set Pi, therefore using proposition 4.2 we have that each point of Pq has two preimages. 
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Figure 6: The union of the pre-critical set of the FLM for (a,e) = (2.75,0.12), from top to 
bottom and left to right we have the union of the k first pre-critical sets, for k = 1,2,3,4. The 
horizontal axis corresponds to 6 and the vertical one to x. 



which will belong to T x [0, 1]. These two preimages form the set P-i. In the top right part 
of the figure we have added the set P-2 to the previous ones. The lower component of P-i is 
below Pi (so it has two preimages), which form part of the set P-2- The upper component of 
P-i intersects Pi. Then only some part of it has preimage in T x [0, 1]. Indeed, the preimages 
form two symmetric arches around Pq. Note that the preimage of Pi is Pq, then the preimages 
of the arch of P_i connecting two points of Pq are two arches connecting two points of Pq. 

When further preimages are considered, higher pre-critical sets are obtained. The number of 
components and the shape of these components depend on the relative position of the previous 
pre-critical set. For example in the bottom part of the figure [g] we have added the sets P_3 (left) 
and P_3 U P_4 (right). We can observe that the two arches of P_2 around Pq give place to two 
pairs of arches of P-s, one around each component of P-i. Some parts of the two arches of P_3 
lay below Pi, then when we consider their preimage they give place to new arches of P_4 around 

PoUP_2. 

We are considering these pre-critical sets because they suppose an obstruction to the reducibility 
of the invariant curve x^^i;. Assume that t/q is an invariant curve by the map. Since its graph is 
invariant by the map we have that if it intersects P-k for some k it will eventually intersect the 
set Pq. Therefore given an invariant curve i/q, a necessary condition for its reducibility is that 

P_fc n Graph(yo) = VA; > 0, 

where Graph(2/o) = {{0, x)eTxR\x = yoiO)}. 
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For example, the FLM for the values a and e of the figure [6] we can see that there is a component 
of P-4 which connects Pq with Pi. Then we have that it can not exists an invariant reducible 
curve above Pq. 

We want to formalize this discussion on the obstruction to the reducibility to give additional 
restrictions on the reducibility of the set Xa,e the invariant set of the FLM introduced at the 
beginning of this section. To do that we need to analyze with some more detail the pre-critical 



set. We have the following improvement of the proposition 4.2 



Proposition 4.4. Consider the FLM with a{l + |e|) < 4 and the post-critical set Pi 



{{9,x) G T X 



-(1 + ecos(27r(^ — w)))}. Consider also Si the set of points of the com- 



pact cylinder T x [0, 1] which are below or in Pi, in other words 
Si = i^{9,x) G T X [0,1] 



a 



X < -(1 + ecos(27r(6' 



■a;)))}. 



Then we have that any point (9, x) £ Si \ Pi has exactly two preimages which are given by 
H^{9,x) and H^{9,x), where 



H± : Si 



X [0, 1] 
9 -co 



\ 2 



=F 



a(l + ecos(27r(6l-w))) J 



Moreover we have that (respectively H-) maps homeomorphically the set Si to the set T x 
[0, 1/2] (and respectively Tx [1/2, 1]). Finally we have that the map preserves the orientation, 
in the sense that, it is monotone with respect to x. On the flip side the map H- reverses 
orientation, i. e. it swaps relative positions (with respect to the x- coordinate). 



This result can be used to describe the pre-critical set with some more detail. By construction 



we have that when a > 2 the constraint ( 10 ) is satisfied, consequently the set Pi is strictly above 



Pq. We can apply the proposition 4.4 to obtain that the set P_i is composed by the union of 



two different sets, i. e. P_i = P(+) U P(— ), where 

P(+) = H+[Co) and P(-) = i7_(Co). 
Moreover we have that P{+) is below Pq and P(— ) is above it. 

Now we can consider further preimages. Since P{+) is below Pq we have that it belongs to 5i 
therefore its preimages are defined. Let us denote by P{++) = H-^{P{+)) = Hj^ o H^{Pq) and 
P(+-) = H.{P{+)) = H. o i/+(Po). 

On the other hand when, we consider the preimages of P(— ), we can have different relative 
positions between P(— ) and Pi depending on the parameters. It might happen that the curve 
P(— ) is completely below Pi, completely above it or that they intersect. In the case that P(— ) 
has points above Pi we have that H± is not defined in these points, therefore the set H±{P{—)) 
is not well defined. This can be fixed if we formally extend H± to T x [0, 1] as 



//,(^,x) = |^(^'-) if(^,x)E5i 
[ (/) otherwise. 



With this extension we can consider the sets P(— it) = H±(P{—)) without problems of definition. 
If the set P(— ) is completely above Pi we will have that P(— it) = 0, and if the set P(— ) is 
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Figure 7: An example for the symbolic codification of the pre-critical sets. We have the union 
of the pre-critical set Pq to P4, and the post-critical set Pi of the FLM for (a, e) = (2.75, 0.12). 
We have indicated the symbolic codes of some of the components of the pre-critical sets. The 
horizontal axis corresponds to 6 and the vertical one to x. 



partially above Pi we have that P(— it) is the preimage of the points below Pi. With this 
notation we have that 

P_2 = P(++) U P(+-) U P(-+) U P( ). 

This symbolic codification of the different components of P_2 has a straight forward generaliza- 
tion to the sets P_fc. 

Proposition 4.5. Let {+, denote the Cartesian product k times of the set of two elements 
{+, —}. Given s = (si, . . . , Sk) G {-|-, , let us define the set P(s) as 

P{s)=H,^oH,^_^o...oHs,{Po). 

Then we have, 

p_fc = u ^(^)- 

In figure [7] we have the four first pre-critical sets of the FLM for the parameter values (a,e) = 
(2.75,0.12). The symbolic codification of some of the components of the pre-critical set have 



been indicated in the picture. The different sequences can be deduced using proposition 4.4 



Now we will see which components of the sets P^k can pose an obstruction to the reducibility 
of the invariant curve Xa,e- 

Consider the pre-critical sets of the FLM for a > 2 and suppose that the first constraint of 



reducibility (10) is satisfied. Then we have that P(— ) is a closed curve above Pq. To construct 
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Figure 8: The different constrains on tlie region of reducibility of tlie invariant curve Xa^e together 



with the bifurcation diagram of section 3.2 The curve (a) corresponds to the constraint (10), 
the curves (b), (c), (d) are the constrains associated to the pre- critical set P((-)^), P{{-Y) and 
)^^) respectively. 



additional constrains on the reducibility let us assume that P{—) and Pi intersect in exactly 
two different points. If we consider the preimages of P{—) we have that P{ — h) is connected 

by an arch to Pq and below it. On the other hand P{ ) is also an arch connected to Pq but 

above it. Now it can happen that P( ) intersects the post-critical set Pi. If such is the case 

then we would have that one of the component of P_2 connects Pq with Pi . This would imply 
that no invariant reducible curve can exists between Pq and Pi. 

To avoid the situation described above, one might require the set P( ) (when defined) to 

be below the set Pi. This condition gives place to an additional constraint in the parameter 
space of the FLM. In figure [8] we have this constraint on the reducibility together with the first 



constraint and the bifurcation diagram of section 3.2 



When the set P( ) exists but it stays below Pi we can consider further preimages of the 

set. We have that P(— ± +) = P'_|_(P(— zt)) will be below Pq since we have that maps 
5i homeomorphically to T x [0,1/2]. Then this set does not suppose an obstruction to the 
reducibility. On the other hand we have that P(— it) define two arches around Co the critical 
set, then we have that P(— =b — ) = iJ_(P(— zt)) define two arches around P(— )• It can happen 
that these arches are below Pi. If this is the case we can consider again their preimages by 
i/_ and H^. The preimages by will be below Pq then they can be discarded, because they 
will not become an obstruction to the reducibility of Xa,e- The set P( — h) is below Pq, then 

P(- + -) = F_(P(-+)) will be above P(-) and then P(- + ) = H-{P{- + -)) will be 

below H-{P{—)) = P( ). If P( ) does not suppose an obstruction to the reducibility, 

neither does P( — I ). Finally the set P( ) will be an arch above Pq U P( ), then 

this can intersect the set Pi becoming an obstruction to the reducibility. 
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To avoid P{ ) being an obstruction to the reducibility one should require it to be below 

Pi for any parameter value. This will be an additional constraint to the previous ones. In figure 
[Hwe have also added this constraint. 

Finally note that the argument used above can be extended to any order. Assume that we 
have P((— )^") which exists but it is not an obstruction to reducibility, where (— represents 
the 2n times repetition of the symbol — . We will have that -ff+(P((— )^")) can be discarded 
for being below Pq, and H+ o can be discarded for being below the union of all 

the set P{{—)'^^) for k < n, then the only set which can suppose an additional restriction is 

o i7_(P((-)2")) = p((_)2(n+l)). 

Note that the higher order conditions does not necessarily suppose an improvement of the 
previous ones. For example in figure [8] we have considered the additional constraints, but not 
until P{{—)^'^) we have had an improvement to the constraint given by P((— )^). 

4.3 Remarks on the constraints 

We have assumed that the sets P((— )'^) for k odd intersect the set Pi in exactly 2 points, but we 
actually have that this assumption can be omitted. If the curves intersect in an even number of 
points we have that the discussion above is still valid but considering the different arches at the 
same time. When they intersect in an odd number of points there is at least one point where 
both curves are tangent. Then the preimage of the intersection is a single point of Pq, therefore 
it can be omitted. On the other hand, we have only taken into account the case where an arch of 
P((— )2^) intersects the set Pi as a possible obstruction to the reducibility of the curve Xa,e ■ It 
might also happen that an arch of the set P{{—)'^^) intersects another arch of the set P{{—)'^^~^) 
for some fc > r > 0. This can produce also an obstruction to the reducibility of the invariant 
curve Xa,e as well. But this case can be omitted because, if it occurs, then we can consider the 
(2r + 2)-th image of both sets by the FLM an we wiU have that p((-)2(fc-»-i)) n Pi / 0. 

Some other properties can be deduced for the pre-critical sets. For example, if P((— )^°) = 
for some ko, then we have that P{{—)^) = for any k > k^. On the other hand we have that if 
P((— )2'^") n Pi 7^ for some ko then P{{—)^) ^ for any k > k^. Another interesting property 
is the fact that when the set P{{—)'^^) does not intersect Pi for any fc > then we have that 
between them there exist an invariant compact subset in T x [0, 1]. The set is delimited on the 
top by the union of the set Pi and the sets P{{—)'^^~^^) for any k > 0, and below by the unions of 
the set Pq and the sets P{{—)'^^') for any r > 0. Note that any curve inside these regions will be 
reducible. Indeed, when this invariant subset exists we will say that the FLM has an invariant 
set of reducibility. 

In the previous subsection the constrains on the reducibility have been given through geometric 
conditions on the pre-critical set. The numerical computation of the constrains has been done 
as follows. We have fixed a value for a, for instance a = 3.2. Then we have localized a value 
of e such that P{{—)'^^) exists. We know that for e = 1 — 1/a this is satisfied, but for a better 
numerical stability of the computation it has been convenient to consider lower values of e. 
Once the parameters have been fixed we have looked for the points P{{—)'^^) H Pi, then these 
points have been continued in e until a tangency between the curves has been obtained. Then 
the tangency has been continued in the parameter space as a constraint on the reducibility. A 
priori, the constraints obtained do not have to be better than the previous ones, therefore only 
the improvements of the constraint have been kept. 
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Figure 9: Some remarkable sets of the FLM for (a, e) = (3.2, 0.138). In a solid line we have the 
union of the sets Pq, P-i, P-2, P-3 and all the sets of the type P{{—)^) for any A; > (recall 
that we admit P{{—)^) = 0. In a dashed line we have the post-critical set Pi and in a dotted 
line the attracting curve of the map. The horizontal axis corresponds to 9 and the vertical one 
to X. 



Analyzing the results on the figure |8l one can observe that for values of a close to 2, the first 



constraint (10) coincides with the boundary of the parameter space where the attractor is a 
reducible curve. When the parameter a gets further from 2 the constraint is replaced by the 
constraint given by the set P((— )^), later on the constraint given by the set P((— )^) becomes 
optimal, and even further this one is replaced by the constraint associated to the set P((— )^^). In 
other words it seems that the minimum of all the possible constrains gives the optimal constrain 
on the reducibility of Xa,e- We think this can be due to the fact that whenever there exists an 
invariant set of reducibility, then the invariant curve Xa,e remains in the set. We have only found 
numerical evidences for this fact. In figure [9] we have plot the attractor of the FLM together with 
pre-critical sets of the map. We can observe that, even though the invariant set of reducibility 
is very thin, the invariant curve lays completely in its interior. 

Then the different constraints on the reducibility would be optimal because they give the region 
(in the parameter space) of existence of the invariant set of reducibility. For parameters close 
to 2 the invariant set of reducibility is delimited only by Pq and Pi. When the parameter a 
gets further from 2, at some point in the parameter space the invariant set of reducibility gets 
delimited by Pq U P((— )^) and Pi U P(— ). Then, from that point on, the constraint which 
determines the existence of the set is the tangency between P((— )^) and Pi, what makes the 
minimum of both constraints the optimal one for the reducibility. When the parameter a gets 
even further from 2 the existence of the invariant set of reducibility is determined by higher 
order pre-critical sets, and the one which determines the existence of the set gives the optimal 
constrain. 

A final remark is that the constrains obtained do not depend on the stability of the invariant 
curve. In other words, they are also valid after the period doubling of the attracting set. In 
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Figure 10: In the top left we have an hivariant repeUing curve x{9) of the FLM ^ for a fixed 
at 3.5 and £ just before losing its reducibility (which happens for e ~ 0.154086). In the top 
right and bottom left we have the first {x'{8)) and the second derivative {x"{9)) of the same 
invariant curve. In these pictures the horizontal axis corresponds to 0. In the bottom right we 
have (in the vertical axis on a logarithmic scale) the supremum |Xg(^)| (solid line) and |x"(0)| 
(dashed line) as a graph of the parameter e (horizontal axis), where Xe is the continuation of 
the invariant curve xq{9) = 1 — 1/a for a fixed at 3.5. 



the unstable case, it is known that the IFT can not be applied to prove the persistence of 
an invariant curve after the reducibility is lost (see section 3.5 of |20]). Indeed we have done 
numerical computations which suggest that when the invariant curve loses its reducibility it 
suffers a "fractalization" process. Concretely, we have continued numerically with respect to 
£ the curve Xa^e for a = 3.5. To compute the invariant curve, we have approximated it by 



its truncated Fourier series as described in section 3.2 Looking at the bifurcation diagram in 



figure |4] we have that this lays completely to the right of the period doubling bifurcation. In 
figure [To] we have displayed the invariant curve and its derivative for the last parameter e of 
the continuation. We can see that the curve looks quite sharp, due to big variations of its 

derivative. In the same figure we have sup ^{9)\ and sup \ x^ ^(0)1 against £. We can observe 

6»eT ' 6»eT 

that that sup jx'^ ^(0)| seems to grow unbounded when the parameters of the map get closer 

to the boundary of reducibility. Note that the process of destruction is the same as the one 
obtained when moving the parameters through the period doubling bifurcation. If an unstable 
curve is really destroyed when it loses its reducibility, then we would have that the constraints 
on the reducibility actually enclose the region of existence of the curve Xa,e in the unstable case. 
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5 Period doubling and reducibility 



In section 3.2 we have studied numerically the bifurcation of the FLM. Concretely we have 
some bifurcation diagrams of the map in figure |4] It can be observed that from each parameter 
value (a,e) = (/i,0) it is born a curve of period doubling bifurcation, which is confined in a 
reducibility zone. Moreover, the period doubling bifurcation curve collides (in a tangent way) 
with the boundaries of the reducibility, giving place to a codimension two bifurcation, namely 
"period doubling - reducibility loss" bifurcation. In the first part of this section we will assume 
that we have an invariant curve and we give a result which says that this tangent collision should 
be expected. In the second part we will give a model map for the reducibility regions observed 
in the bifurcations diagram of the FLM. 



5.1 Interaction between period doubling and reducibility loss bifurcations 

Theorem 3.3 of |20j describes the behavior of the Lyapunov exponent A = A(/x) of a one para- 
metric family of maps, with respect to its parameter fi close to a value /io where the reducibility 
is lost. Now we will present a result which have some resemblances with the cited one, but in the 
following result we consider a two parametric family of maps, and then we study the interaction 
of the period doubling and reducibility loss bifurcations in the parameter plane. 

Theorem 5.1. Consider a one dimensional linear skew product 

= e + u, 1 

X = a{6,fj,,X)x, J 

which depend (C°° -) smoothly on 9 and (also (C^ -) smoothly) on two parameters {fx, X) £ U with 
U cM? an open set. 

Suppose that there exists a regular curve in the parameter space, parameterized by {iJ-oit), Xo{t)) 
for t G [0, eo); such that the skew product undergoes a loss of reducibility when we cross transver- 
sally that curve. Let A(/i, A) denote the Lyapunov exponent of the invariant curve. Suppose also 
that A(/io(0), Ao(0)) = and |A(/xo(t)), Ao(t))t=o / 0. 

Then there exists a curve in the reducibility zone of the parameter space, parameterized by 
(/zi(t), Ai(t)) fort G [0, ei), with a quadratic tangency with the previous curve and such that 

A(/ii(t),Ai(t)) = 0/or anyt£ [0,ei). 

Proof. Consider a change of variables to new parameters (/i. A) such that the curve (yUo(i)) '^o(*)) 
goes to the positive // semiaxis. Then the reducibility loss occurs when A crosses transversally 
zero, for any value for fi G [0, e). We can apply the normal form near a reducibility loss (Lemma 
3.5 of |20j ) for any /2 when A = 0. Then we have that there exist a change of variables such that 
the system (llip can be transformed to the following one 



= O + u, 

y = /i(z.(A))(KA) + cos(27r(0-0o(A))))y, 



(12) 



where i^(A) is a smooth function of the parameter A, with 

1/(0) = 1, (13) 
^KO) / 0, (14) 
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and h{i') is a smooth function which never vanishes and 9q{X) a smooth function to T. 



Note that when applying Lemma 3.5 of |2D] above, the original map depends smoothly on the 
parameter /Z. Then the change of variables will also depend on it and so will do the function 
u, h and ^o- Then we can think on these functions depending on both parameters and write 
u = f^(A) = A), h{i^) = hfi{vfi(\)) = h{fl, A)) and 6o{\) = 9o{fl, A). 



The condition (14) must be satisfied for any /j, G [0,e), then ^z^(//,0) can not vanish. Changing 



/i by —fi if necessary, we can suppose that 0) is positive. Therefore conditions (13) and 



(14) can be rewritten as 



1 for any // G [0, e), 



dX 



iy{fj.,0) > for any /i G [0, e). 



(15) 
(16) 



Note that we are assuming that we have reducibility for A > and non-reducibility for A < 0. 



Using the formula ( 25 ) we can compute the Lyapunov exponent of the system, 

"1 



A(/i,A) 



In \hifl, u{fl, A))(z^(/2, A) + cos(27r(6' - Ooifl, A))))|d6' 







/ ln\h{fl,u{jl,X))\de+ I + cos{2tt{9 - eQ{li,X)))\de 

Jo Jo 



In 



In 



h{^i,v{n,\)) 



if A < 0, 

+ arccosh |z^(/i, A)| if A > 0. 



(17) 



By hypothesis we have that ^A{fj,Q{t)), Ao(i))t=o 7^ 0, performing the change of variables to the 
original parameters, we have that 

d - d 

^A(/2,A)(^_3,)=(o,o) = j^HMt)),Ht))t=o + 0. 



If we differentiate equation (|17|) with respect to \x with the constraint A = we have 
d 



|. {h{ji, 0)) _ d-^h{jL, 0) + dMy{T^. 0)) • d-^v{ji, 0) 



Differentiating in (15) we have (?^i^(;U, 0) = 0, then from the last two equations it follows that 

a^/i(o, 1/(0,0)) 



/o. 



(18) 



Note that for A > we have that 



A(//, A) = if, and only if, |i/(;U, A)| = cosh ( In 



\h{fi,i^{fj.,X))\ 
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As i'{fi,0) = 1, by continuity we have \i'{p,,\)\ = i>{p,,X) for any A small enough. But then, 
recall the h{fi, A)) is a smooth function which never vanishes (in a neighborhood of A = 0). 
Then, if we consider the function 

\h{ll,i'{jl,X))\\ ^^g^ 



A) := A) — cosh ^In 

we have that it is a smooth function on the parameters (/i. A) for G [0, e) and A small enough. 
Moreover from (17) we have that (when A < 0) A(/i, A) = if, and only if, p(//. A) = 0. 

To prove the theorem we first will apply the IFT to the function p{fl, A), then we will check that 
the curve that we obtain is tangent to the set A = and finally that it belongs to the upper 
semi-plane A > 0. 



From equation (17) and the hypothesis A(0,0) = it follows that 



(20) 



Therefore using (15) it follows that p{0,0) = 0. 



To apply the IFT we need to compute first the derivative of p with respect to A and check that 
it is different from zero. We have 

d 



dX 



p{p, A) = d)^i'{p,, A) — sinh I In 



\h{p, v{p, X))\ \ duh{p, u{p, A)) • d-xu{p, A) 



hin,iy{iJ,,X)) 



Evaluating at (//,A) = (0,0) and using (16) we have 

d 



dX 



p(o,o) = a^i/(o,o)>o. 



such that A(0) = 
(21) 



Then we have that there exist an interval [0, ei) and a function A : [0, ei) - 
and A{p, X{p)) = for any p G [0, ei). Moreover we have that 

To finish with the proof we only have to check that A'(0) = and A"(0) > 0. Differentiating in 

\h{p, v{p, X))\ \ djih{p, X) + dyh{il, iy{p, A)) • dfii^{p,, A) 



(19) with respect to p we have. 



d_ 
dp 



p{p, A) = dfii'{p, A) — sinh I In 



h{p,u{p,X)) 



(22) 



Using equations (16) and (20) we can evaluate the last equation at {p,X) = (0,0) and we have 



^^p{p,X) = 0, therefore (using (21)) we have X'{p) = 0. 



For the condition A"(0) > 0, differentiating on equation (21) and evaluating at (and using also 
X'{p) = 0) we have that 



A"(0) 



dfi' 



■p(0,0) 



9aP(0,0) ■ 



Finally, from ( 15 ) we have also that i^i/(0, 0) = 0, then differentiating equation (22 ), evaluating 
at {p, A) = (0, 0) and simplifying all the zero terms we have 



dp 



2P(0,0) 



dfih{o,u{o,o)) 

/i(0, 1/(0,0)) 



which is different from zero due to (18), therefore A"(0) > 0. 



□ 
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5.2 A model for the reducibility regions 



In this section we present a two parametric skew product to model the interaction between the 
reducibihty loss and the period doubling bifurcation We will see that this map is helpful for the 
understanding of the reducibility regions observed in the parameter space of the FLM. 

Let us recall first the period doubling (and pitchfork) bifurcation for one dimensional maps. 
Consider the one dimensional map x = x{x'^ — fi), where /i is a parameter in the real line. It 
is well known that this map undergoes pitchfork bifurcation when the parameter crosses the 
value —1 and it undergoes a period doubling bifurcation when n crosses the value 1. 



The model map that we propose is the following, 

(23) 



X = x{x'^ - {fi + X cos{2tt9))) 
where n and A are parameters, and uj is Diophantine. 

Note that when A is equal to zero the system uncouples and we obtain the model for the generic 
unfolding of the period doubling bifurcation when fi crosses 1. Indeed the q.p. forcing has 
been considered in such a way that the set x = 0, namely the trivial invariant set, is always an 
invariant curve of the map (for any parameter values). 

The linear dynamics around this trivial invariant set is 

9 = e + u:, 1 

X = -(/i + Acos(27r0)). J ^ ' 



We are in the C°° framework, therefore by Corollary 1 in [20] we have that the trivial invariant 
set is reducible if, and only if, |A| < In other words, the loss of reducibility bifurcation 
correspond to the lines in the parameter space given by A = it^. 

On the other hand, we can compute A = A(/i, A) the Lyapunov exponent of the trivial invariant 
curve X = 0. Recall that, 

'l„|r + c„s(2,«)M<,= {:;°^^_ I (25) 



and arccosh |r| = In ^|r| + x/r^ — 1 j for |r| > 1. Using this it is easy to check that 

A(^,A) = / ln|^ + Acos(27r6l))|(i6' 
Jo 



In ^ if |A| > ImI 



Concretely, for |A| < \^\ we have that the skew product associated to the trivial invariant set 
X = is reducible to 

9 = 9 + uj, 



H + sign(;u)x//i2 - A2 



X. 
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Figure 11: Different remarkable curves of the map (23) are displayed. In a solid line we have 



the change of stability of the trivial invariant set. In a dashed line the loss of reducibility of the 
attracting set, which is the trivial invariant set or a two periodic invariant curve. Finally we 



have in a dotted line the boundary of applicability of the proposition 5.2 The horizontal axis 
corresponds to the parameter and the vertical one to A. 



Therefore, in the reducible case (|A| < the trivial invariant set changes its stability when 

l/^l = 1 + -r- 



One should expect that the map ( 23 ) has a pitchfork bifurcation when one crosses the parabola 
/i = — 1 — ^ and respectively a period doubling bifurcation when crosses /x = 1 + Unfor- 
tunately this can only be proved in a small region of the reducibility zone. More concretely we 
have the following results. 



Proposition 5.2. The model map (23) has 



Two-periodic (continuous and non-trivial) invariant curves when < |A| < and 1 + ^ < 



i-2|A|. 



Two (continuous and non-trivial) invariant curves when // < — |A| < and + 2|A| < 
H<-1 + ^. 



Proof. We will prove only the first item of the proposition, since the other one is completely 
analogous. 

We are interested on the two periodic invariant curves of ( |23[ ). Let us denote by / the function 
which defines the system, i.e. f{0,x) = x (x^ — (/x + A cos(27r6'))) . Observe that f{9,—x) = 
—f{6,x) for any {9,x) G T x M and any parameters values. It follows that, if : T — )• M is an 



invariant curve of (23), then —u also is. Moreover, given a curve u (different from the trivial 



invariant set u = 0) satisfying 

u{e + w) = -f{9, u{e)), for any OgT, 
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we have that u is a two periodic solution. 



Let us consider the system 



X 



e + oj 

h{9,x), 



(26) 



where h{6, x) 



-f{0,x) 



-X (x^ 



(;U + A cos(27r^))) . It is clear that an invariant curve 



(different from x{9) = 0) for the system above is a two periodic solution of (23). We will apply 



theorem 4.2 of [15j to this map in order to prove the existence of periodic solutions of (23). Let 
us check the conditions of the theorem. We have that 



d_ 
dx 



h{9, x) = — 3x^ + ^ + A cos(27r6 



When < |A| < /i we have that ^h{9,x) > for any 
Concretely this condition is satisfied when < ^— 



, x) such that x^ < A'+-^cos(27re 



On the other hand, we have that 
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2/1(6*, x) 



therefore 



dx 



s^h{e,x) 



Q3 

-6x and -—^h(9,x) 

OX'^ 



-6 



-6x 



lh{e,x) 2 \lh{9,x) 



Note that whenever ^^h{9, x) > then we have that Sxh{9, x) < 0. Finally the only hypothesis 
remaining to be satisfied is to find an interval I such that T x I is invariant by the system ( 26 ) . 



Let us consider the interval / 



M-|A| 



3 



S 



with S an arbitrary small value yet 



to be defined. Note that for any point on the interval the monotonicity condition is satisfied, 
then it is enough to check that the interval satisfies the invariance condition. In order to have 
invariance of the interval we have to check that 



\h{9,x)\ < 



^-|A| 



for any {9, x) G T x /. Since the monotonicity condition is satisfied, it is enough to check that 



/i-|A| 



6 



< 



f,-\X\ 



5. 



Then it follows that it is enough to have 

2^ ^|A| 



+ Acos(27r6i) 



< L 



when /i + 2|A| < |, we have that there exists a value of S sufficiently small such that this is 
satisfied for any 9 £T. 

We are in situation of applying the theorem 4.2 of [15j, we know that u{9) = is always a 
continuous invariant curve, which is contained in the set T x /. Moreover we know its Lyapunov 
exponent explicitly, therefore when this crosses zero, the theorem implies that there exist two 



invariant curves (with negative Lyapunov exponents) of the system (23), which correspond to 



periodic solutions of the system (26). Moreover the invariant curves have negative Lyapunov 
exponent, therefore the periodic invariant curve of the original map is attracting. □ 
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To illustrate this last proposition, in figure 11 we have plotted the curve which constrains the 
validity of proposition 5.2 



Let us assume that in the reducible case (|A| < |/i|) the parabola /i = 1 + ^ corresponds to 
a period doubling bifurcation. Note that this parabola has a tangency at the points (/u, A) = 



(2, ±2) with the boundary of reducibility, as predicted by theorem 5.1 Then these points in the 
parameter space would correspond to the "period doubling - reducibility loss" bifurcation, since 
it is the point where both curves merge. 



Recall that in section 3.2 we have reported how the period doubling bifurcation of the FLM were 



enclosed inside regions of reducibility. In the proposed model (24) for the reducibility regions it 
is easy to justify this behavior. 

Differentiating the equation which defines the map we have that the critical region of the map is 
Pq = {{9, x) G T X /| 3x^ — (/i + A cos(27r0)) = 0}. Following the arguments of sectionjZ] we have 



that a reducible invariant (or periodic) curve of (24) can not intersect the critical set. Note that 



for |A| < In] the set Pq is composed by two closed curves in the cylinder, one in each side of the 
trivial invariant set. Moreover when |A| tends to |^| we have that the components get closer to 
the trivial invariant set x = 0. 



If we assume that the parabola in the parameter space // = 1 + ^ corresponds to a period 
doubling bifurcation, then we have that when the trivial invariant set becomes unstable then 
a period two solution must be created in its neighborhood. But then, when |A| tends to 
we have that the set Cq gets closer and closer to the trivial invariant set, therefore there is no 
room for the period doubled curve to be reducible. This explains why arbitrarily close to the 
"period doubling - reducibility loss" bifurcation parameter (/x, A) = (2, ±2) one can observe a 
reducibility loss of the period two invariant curve. 



In figure 11 there are shown the different bifurcation curves of the map (24). The reducibility of 
the period two invariant curve have been estimated numerically. Let us remark the resemblance 
of the region of reducibility of the attractor with the same regions of the bifurcation diagram of 
the FLM. 



6 Summary and conclusions 

We study the Forced Logistic Map as a toy model for the truncation of the period doubling 
cascade of invariant curves. In section [33] we have done a numerical analysis of the bifurcation 
diagram of the FLM, which is displayed in figure |4j This computation revealed that each 
period doubling bifurcation curve in the parameter space is confined inside a region where the 
attracting invariant curve is reducible. Now we can use studies done in sections |4] and [5] to 



review the analysis of the bifurcation diagram of the FLM done in section 3.2 



In section [4] we have done a study of the critical set, their images and their preimages. We have 
constructed different constrains in the parameter space for the reducibility of the invariant curve 
Xa,e (which is the continuation of the invariant curve Xa,o{0) = 1 — 1/a for e > 0). We have 
also illustrated how the combination of all these constrains seems to be the optimal constrain 
for the reducibility of the invariant curve. In other words, they approximate the boundary of 
reducibility of the attracting set of the map. Using the notation introduced in section 3.2, we 
have that these constraints characterize the curve Cq . In the bifurcation diagram of the figure 
[4] only the properties of the stable set are reflected, but we have that the constrains are still 
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valid after the period doubling. Actually we conjecture that these constrains give the boundary 
of existence of the curve Xa,e when it is unstable. 

In section [5] we have studied the interaction between the reducibility loss and the period doubling 



bifurcation. For the case of linear skew products we have theorem 5.1 which says that generically 
we can expect the period doubling bifurcation curves and the reducibility loss bifurcation curve 
to be tangent. The diagram of the figure |4] has been done in terms of the attracting set of the 
FLM. As the FLM is not a linear skew product, this theorem is not applicable. But it can be 
applied to the linear skew product given by the linearization of the map around the invariant 
curve. In the same section, we have also given a model for the interaction of the period doubling 
bifurcation and the reducibility loss. With this model we have seen that, if the period doubling is 
close to a reducibility loss, then there is an obstruction to the reducibility of the period doubled 
curve. This explains why the curves Cq , Di and meet at the same point (using again the 



notation of section 3.2). Moreover, theorem 5.1 gives a good explanation of why they do it in 



a tangent way. Finally, let us remark that the study done in section [5] does not depend on the 
map considered, therefore it can be extended to the rest of reducibility regions determined by 
C^_i and C~ (and containing Di). 

Part of the study done here will be continued in |17 1 I18 1 [T9]. In |17j we will propose an extension 
of the renormalization the theory for the case of one dimensional quasi-periodic forced maps. 
Using this theory we will be able to prove that the curves of reducibility loss bifurcation 
really exists (for e small enough). In [18] we will use the theory proposed in the previous one 
to study the asymptotic behavior of the reducibility loss bifurcations when the period goes to 
infinity. In the previous two articles several conjectures will be done. In [TH] we will support 
numerically this conjectures. We will give also numerical evidences of the self-renormalizable 
character of the bifurcation diagram of the figure |4] when different values of the rotation number 
u are considered. 



References 

[1] R. A. Adomaitis, I. G. Kevrekidis, and R. de la Llave. A computer-assisted study of global 
dynamic transitions for a noninvertible system. Internal. J. Bifur. Chaos Appl. Sci. Engrg., 
17(4):1305-1321, 2007. 

[2] L. Alseda and S. Costa. On the definition of strange nonchaotic attractor. Fund. Math., 
206:23-39, 2009. 

[3] A. Arneodo, P. H. Coullet, and E. A. Spiegel. Cascade of period doublings of tori. Phys. 
Lett. A, 94(l):l-6, 1983. 

[4] K. Bjerklov. SNA's in the quasi-periodic quadratic family. Comm. Math. Phys., 286(1):137- 
161, 2009. 

[5] H. W. Broer, G. B. Huitema, F. Takens, and B. L. J. Braaksma. Unfoldings and bifurcations 
of quasi-periodic tori. Mem. Amer. Math. Soc., 83(421):viii+175, 1990. 

[6] H. W. Broer, C. Simo, and R. Vitolo. Bifurcations and strange attractors in the Lorenz-84 
climate model with seasonal forcing. Nonlinearity, 15(4):1205-1267, 2002. 

[7] H. W. Broer, C. Simo, and R. Vitolo. The Hopf-saddle-node bifurcation for fixed points 
of 3D-diffeomorphisms: the Arnol'd resonance web. Bull. Belg. Math. Soc. Simon Stevin, 
15(5, Dynamics in perturbations): 769-787, 2008. 



31 



E. Castella and A. Jorba. On the vertical families of two-dimensional tori near the triangular 
points of the bicircular problem. Celestial Mech. Dynam. Astronom., 76(l):35-54, 2000. 

U. Feudel, S. Kuznetsov, and A. Pikovsky. Strange nonchaotic attractors, volume 56 of 
World Scientific Series on Nonlinear Science. Series A: Monographs and Treatises. World 
Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2006. Dynamics between order and 
chaos in quasiperiodically forced systems. 

J. Figueras and A. Haro. Computer assisted proofs of the existence of fiberwise hyperbolic 
invariant tori in skew products over rotations. To appear, 2010. 

V. Franceschini. Bifurcations of tori and phase locking in a dissipative system of differential 
equations. Physica D Nonlinear Phenomena, 6:285-304, April 1983. 

A. Haro and C. Simo. To be or not to be a sna: That is the question. Available at 
|http : / / www . maia . ub.es/ dsg/2005/ ind ex . html , 2005. 

J. F. Heagy and S. M. Hammel. The birth of strange nonchaotic attractors. Phys. D, 
70(1-2):140-153, 1994. 

M. W. Hirsch, C. C. Pugh, and M. Shub. Invariant manifolds. Bull. Amer. Math. Soc, 
76:1015-1019, 1970. 

T. H. Jaeger. Quasiperiodically forced interval maps with negative schwarzian derivative. 
Nonlinearity, 16(4):1239-1255, 2003. 

A. Jorba. Numerical computation of the normal behaviour of invariant curves of n- 
dimensional maps. Nonlinearity, 14(5):943-976, 2001. 

A. Jorba, P. Rabassa, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps I. Existence of reducibily loss bifurcations. In 
preparation, 2011. 

A. Jorba, P. Rabassa, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps II. Asymptotic behavior of reducibility loss bifur- 
cations. In preparation, 2011. 

A. Jorba, P. Rabassa, and J.C. Tatjer. Towards a renormalization theory for quasi- 
periodically forced one dimensional maps III. Numerical evidences. In preparation, 2011. 

A. Jorba and J. C. Tatjer. A mechanism for the fractalization of invariant curves in quasi- 
periodically forced 1-D maps. Discrete Contin. Dyn. Syst. Ser. B, 10(2-3) :537-567, 2008. 

K. Kaneko. Doubling of torus. Progr. Theoret. Phys., 69(6):1806-1810, 1983. 

S. Kuznetsov, U. Feudel, and A. Pikovsky. Renormalization group of scaling at the torus- 
doubling terminal point. Phys. Rev. E (3), 57(2, part A):1585-1590, 1998. 

A. S. Pikovsky and U. Feudel. Characterizing strange nonchaotic attractors. Chaos, 
5(l):253-260, 1995. 

A. Prasad, V. Mehra, and R. Ramaskrishna. Intermittency route to strange nonchaotic 
attractors. Phys. Rev. Let, 79(21):4127-4130, 1997. 

A. Prasad, V. Mehra, and R. Ramaskrishna. Strange nonchaotic attractors in the quasiperi- 
odically forced logistic map. Phys. Rev. E, 57(2): 1576-1584, 1998. 



32 



[26] P. Rabassa. Contribution to the study of perturbations of low dimensional maps. PhD 
thesis, Universitat de Barcelona, 2010. 



[27] C. Simo. On the analytical and numerical approximation of invariant manifolds. In Les 
Methodes Modernes de la Mecdnique Celeste (Course given at Goutelas, France, 1989), D. 
Benest and C. Froeschle (eds.), pages 285-329. Editions Prontieres, Paris, 1990. Available 



at http : / /www . maia . ub . es/dsg/2004/ index . html[ 



[28] L. van Veen. The quasi-periodic doubling cascade in the transition to weak turbulence. 
Phys. D, 210(3-4):249-261, 2005. 



33 



